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A study of the effect of unscreened Coulomb interactions on the quantum spin Hall (QSH) phase 
of finite- width neutral zigzag graphene ribbons is presented. By solving a tight-binding Hamiltonian 
that includes the intrinsic spin-orbit interaction (I-SO), exact expressions for band-structures and 
edge-states wavefunctions are obtained. These analytic results, supported by tight-binding calcula- 
tions, show that chiral spin-filtered edge states are composed of localized and damped oscillatory 
wavefunctions, reminiscent of the ones obtained in armchair ribbons. The addition of long-range 
Coulomb interactions opens a gap in the charge sector with a gapless spin sector. In contrast to 
armchair terminations, the charge-gap vanishes exponentially with the ribbon width and its am- 
plitude and decay-length are strongly dependent on the I-SO coupling. Comparison with reported 
ab-initio calculations are presented. 
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Graphene, the single-layer carbon material isolated for 
the first time in 20040, i provides a unique window 
to explore phenomena usually associated with QED as 
well as quantum Hall physics 0, 0]. Its peculiar band- 
structure, with linear spectra at specific points in mo- 
mentum space, allows for a description of low-energy 
properties in terms of a Dirac Hamiltonian for massless 
electrons. Besides the interest on understanding its fun- 
damental properties, graphene has received increasing at- 
tention due to its potential technological applications. In 
particular, much effort is focused on obtaining graphene 
nanoribbons (GNR) with tailored properties [5|, |6|, 13 that 
could be key elements in future electronics and spin- 
transport applications 

Ribbons present an exceptional opportunity to study 
confinement effects and enhanced Coulomb interactions 
on graphene properties. Moreover, they arc a natural 
test-bench to investigate the physics of the quantum spin 
Hall phase (QSH) originated by intrinsic-spin orbit (I- 
SO) interactions [9(. This phase has been studied in the 
low-energy limit (Dirac Hamiltonian) and was found to 
be characterized by a bulk insulating behavior and chiral 
spin-filtered edge states. Tight-binding and perturbative 
renormalization group calculations 0, lloL 11 1 , predict the 
phase to be stable against disorder and screened Coulomb 
interactions along the edges. These arguments also sug- 
gest a transition to an insulating phase for unscreened 
interactions. However, at present no detailed study ex- 
ists on its properties and the fate of the spin- filtered edge 
states when the full graphene band-structure is consid- 
ered and inter-edge unscreened Coulomb interactions are 
included. 

In this letter we present the first study of the QSH 
phase in zigzag GNR (ZGNR) that includes the com- 
bined effects of finite-system size and electron-electron 




FIG. 1: A zigzag graphene nanoribbon (ZGNR) (left) of width 
W and its corresponding Brillouin zone (right). Dotted lines 
indicate positions where boundary conditions are applied. 



interactions. The analytic solution of a tight-binding 
Hamiltonian with I-SO interactions provides exact ex- 
pressions for the band-structure and wavefunctions that 
characterize the QSH phase giving insight into the nature 
of the spin-filtered edge states. To investigate the effects 
of unscreened Coulomb interactions we exploit the linear 
nature of the spectrum and use bosonization methods to 
analyze the Hamiltonian. Our results indicate that a gap 
opens in the charge sector while the spin sector remains 
gapless. Within the harmonic approximation, we pro- 
vide an expression for the gap in terms of the parameters 
of the model. The analysis of the data suggests an ex- 
ponential vanishing gap with increasing ribbon's width, 
with an amplitude and decay-length dependent on the 
I-SO coupling constant. These findings are in contrast 
with results for armchair ribbons [12] and provide an in- 
terpretation for reported DFT calculations [3]. 

The Hamiltonian: In the tight-binding description, 
the Hamiltonian Hq = J2<ij>[^ c l c j + ^- c ] (*>J label po- 
sition and spin degrees of freedom); is written in a spin- 
dependent basis, the pseudo-spin basis \l/^ s = (u s A ;u s B ), 
where the left (right) component represents the ampli- 
tude of the wavefunction at a lattice A(B) site, and 



s = ± labels the spin. The intrinsic spin-orbit inter- 
action involves spin-dependent second-neighbor hopping 
terms: Hi-so ~ it'(dik x d>kj)c\s z Cj, with as firsl 
neighbor lattice vectors connecting electrons at positions 
i and fe, cj the electron creation operator at site i anc 
s 2 the spin[^, [f2|- This term satisfies all the symme- 
tries of the graphene lattice 1J] and plays a crucial role 
in magnetic properties of finite-width ribbons. Reportec 
values for t' from ab-initio calculations range from 1.2K 
to lOmK In reciprocal space the Hamiltonian reads 
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where ~/(k x ,k y ) = 2t'(smk x a — 2 sin 2|2 cos k y b), 
4>{k x ,k y ) = t( c ^ 2fc / 3 + 2cos^e-^ fc / 3 ) ! 4>{k x ,k y ) = 
4>(k x , —k y ), and 6 = ay/3/2. The corresponding eigenval- 
ues and eigenstates are E = ±s = ±yj (j>(/> + j 2 and: 
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where r]±(k x ,k y ) = (±e— s~f)/<f> = (j>/ (±e + s^) . N is the 
normalization factor. A gap opens at the Dirac points 
and the system becomes a bulk- insulator Q. 

To study confined geometries we consider a ribbon of 
length L along the x— direction, with finite width W 
along y (see Fig.fTT])), and impose hard- wall boundary 
conditions on the wavefunctions of Eq.(TT]). In addition 
to the set of boundary conditions used in zigzag ribbons: 
(ip A (0) = <p B (W- 6/3) = 0)0, US BE!, the I-SO in- 
teraction introduces new ones: ipa{W) = <^b(— 6/3) 
i i . The more general wavefunction is a superposition 
of four degenerate states of the form of Eq.([2]), with 
k y = (±fci; ±fc 2 ). To guarantee their degeneracy, the pa- 
rameters must satisfy the condition: cosfci6 + cosfc 2 6 = 
2c(l - t 2 ) with r 2 = I/(8t' 2 sin 2 %^). A wavefunction 
for a right mover with spin up is given by: 



* +T = N e ik * x 
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with rjj^ = rj{k y = k/)-\ and N is the normalization fac- 
tor. The new dispersion relation and the coefficients ap- 
pearing in Eq.Q are obtained from the condition: 
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FIG. 2: Band structure of ZGR with W — 46 in the presence 
(left panel) and absence (right panel) of I-SO interactions with 
t' — O.lt. Top inset shows the change in the bulk-bands. Full 
lines are numerical tight-binding calculations. 



with t' = O.lt, in the range ir < k. x < 2ir. The quasi- 
degeneracy of the lower band is lifted for n < k^^ 47r/3 
by the I-SO interaction as predicted [H [H E, EH • The 
spectrum remains gapless at k x = tt, with linear dis- 
persions and right- and left-mover states for each spin 
component The nature of the edge states depends 
on k x : in region / (near k x = tt), the wavenumbers 
ki/2 are complex, rendering damped oscillatory wave- 
functions; in region 77 (extending to Ak x ~ 1/W be- 
fore the Dirac points) however, the wavefunction becomes 
fully localized at the edges. The predicted 'spin-filtered 
edge states' are formed by a combination of states from 
these two regions. As k x is increased further, fc 2 - * 
while k\ is purely imaginary, resulting in oscillatory wave- 
functions in region III . A similar analysis reveals that 
bulk conduction bands are made from a combination 
of localized and oscillatory wavefunctions, except near 
k x = tt where wavefunctions are fully localized (see inset 
in Fig.Q). Notice that a simplified expression for the 
lower-band dispersion can be found when i'/i S> b/W: 

e = ±6t' sm(k x a)t 



1 1 + 1 6i' 2 sin 2 (k x a/2) which near k x = 
tt can be approximated by e w 6t'k x a with the velocity 
of right (left) movers given by hv = ±6t'a. These re- 
sults hold for ribbon widths W = 2Mb with M integer. 
For W = (2M + 1)6 the hard-wall boundary condition 
breaks one of the sublattices discrete translation invari- 
ance symmetry along the y— direction, which opens a gap 
A w i(2*7*) w ' _1 at k x = tt and t'/t < 1. 

The coefficients of the wavefunction (Eq. (J3j)) are given 
by a(fci, fc 2 ) = /V4 e -ifeiW/2 where 

= (f}[ -y 2 )(v'i -r7 2 )(e + 7i)(sinfc 2 6sinfc 2 FF) 
K-?72)W-^)(e + 72)(sinfci6sinfc 1 VF) 1 ' 



1 . = 1 ( ky = and k± = ( fci ± witn ^ = % . e -ifei6/3 ) b = a (-k 1 ,k 2 ),c = a(fc 2 ,fci) and 

d = a(— fc 2 , fei). Representing the resulting spin-up right- 



where q 

fc 2 )/2. The degeneracy condition plus Eq.((3]) uniquely 
determine them as functions of t, t' and W. Fig.@ 
shows the band-structure of a ribbon of W = 46 f» lnm 
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{tp A {y),(pB(W -y)) the 
corresponding spin-up left-mover wavefunction takes the 
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FIG. 3: Probability distribution for a spin-up right-mover for 
values of k x in region I: n < k x a <C n + t' /t; II: n + t' /t <C 
k x a < 4tt/3 - b/W, and III: k x a > 4tt/3 - b/W. Curves for 
spin-down right-movers are obtained reflecting through the 
y = W/2 plane (y/b ->• (W - y)/b, W = 46). 

form *1 T = *[ T = e - ik **(-ip B (W - y),<p A (W - y)) 
and those for spin-down right/lcft-movcrs are obtained 
by taking t' — > — t'. Fig. ([3]) shows the probability dis- 
tribution for a spin-up right (spin-down left) mover as 
a function of the position y across the ribbon for differ- 
ent values of k x . The predicted localized [t| wavefunction 
along the edge changes from exponential to damped os- 
cillatory as k x moves towards k x = n. This oscillatory 
behavior for k x ~ 7r is reminiscent of the one obtained 
in armchair ribbons (AGNR) with I-SO interactions sug- 
gesting an effect independent of particular edge termina- 
tion [121 ]. Also, notice that there is no spin-accumulation 
along the edges neither net magnetization in equilibrium. 

Coulomb interactions: For neutral isolated 
graphene ribbons, screening of the Coulomb potential 
is suppressed due to low electron densities. Starting 
from the un-screened Coulomb potential U pp i (Ax, Ay) = 

e 2 /k^/oq + Ax 2 + Ay 2 p , with p = + (p = — ) labeling 

sub-lattice A (B) and Ay = y — y' + 5 p .- p :a/ V3 [2(3|; the 
interacting Hamiltonian is reduced to 

H=\Y,J d 2 xd 2 yU ± pl iS (x,y)p ± , s ,(x>,y') (6) 

where p±, s = pAs ± Pb s with p ps = c ps c ps and U ± = 
Uaa^Uab- The linear-dispersion of the lower-band near 
k x ~ 7r, makes possible to use the bosonization technique 
to study intra-band Coulomb scattering on spin-filtered 
edge states. Near k x — n, the lower (edge) bands are 
separated by (AE ~ t ~ 2AeV) from the upper bands. 
This gap reduces near the Dirac points to AE ~ l/W but 
does not vanish and, at low enough temperatures, inter- 
band transitions can be neglected. We introduce a spin- 
dependent Bogoliubov transformation from the original 
site-operators to a right/lcft-movers basis: 

CAi(k x ,y) = i PA(k x ,y)i'B^ - fB(k x ,w - y)i>Lt 

CBi(k x ,y) = tpB(k x ,y)ipRi + ifiA{k x ,W - y)i^L\ (7) 
where Vh/lt = ipR/L\(k x ) (for spin-down: (f A /B(y) is 



replaced hy ip B / A (W — y)) . Unlike previous works 0, 
|21| this model includes scattering processes between the 
edges represented by a marginal operator term in the 
Hamiltonian. As a consequence, the non-local irrelevant 
terms are neglected and the resulting Hamiltonian is: 

n = ^Y, [j^(d^c /s ) 2 + K c/s (d x Q c/s ) 2 ] 

c,s c / s 

+ — [cos V87r$ s - cos(%/8tt$ c + Ak F x)] (8) 
na 

where K c/s = yj(l + Q c/s /hv)/(l + R c/s /hv); v c/s = 
v^/(l + R c f s /hv)(l + Q c / s /fiv), and kp is the Fermi mo- 
mentum. The parameters in these expressions are: 

R * = I § E UZ(Ay)B-(y)B"(y') 

J u=± 

R o = J^Y.K{Ay)Fl{y)Fl(y')-R s 
Q> = I E K{Ay)F-(y)F^y') (9) 
with Q c = 0, m = R s /2a, Ay — (y — y') and 

K(y) = l;[<fA{y)(pA(y) + uip B (y)ip B (y) 

+vip B {W - y)(p B (W -y) + uvip A {W - y)ipA(W - y)\ 

B u {y) = ^ A (y)VB(W -y)- u VB {y)vA(W - y) (10) 

In the expressions above, Uq (zero Fourier component) is 
the unscreened part of the Coulomb interaction. Eq.© 
contains forward scattering terms (F™) that preserve chi- 
rality and modify charge and spin- velocities v c / s and Lut- 
tingcr parameters K c / S . Backward scattering terms (B u ) 
originate from the mixing of same spin right- and left- 
movers located at opposite edges and produce the mass 
term to. Using standard approximations [22l |. we obtain: 

A= tvv 1 Zam\ 
a \ nv c J 

Fig((3]) shows the dependence of the gap as a function of 
W for different values of t' as well as the Luttinger pa- 
rameters K c / S indicators of charge and magnetic orders. 
The data was obtained taking k = 2.45 (appropriate for 
graphene over a S1O2 substrate) and ao ~ a. An analy- 
sis of the width-dependence reveals an exponentially van- 
ishing gap A ~ Ac~ BW (see inset in Fig. (J4J)) with gap 
amplitude A and decay-length 1/B decreasing with in- 
creasing t' . These results can be understood in terms of 
the renormalization of the Coulomb interaction, as mea- 
sured by the dimensionless parameter g = e 2 / nhv, with 
v determined by t' . For t' = our results agree with pre- 
vious workfljj. As t' is increased, same spin right- and 
left-movers wavefunction overlaps are suppressed. The 
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FIG. 4: Left pannels: Charge gap A (in units of eV) as a 
function of ribbon width W for a zigzag ribbon of length 
L ~ 25/xm (upper). The armchair (lower) ribbon gap is plot- 
ted for comparison (length L ~ 400pm). Inset in ZGNR 
pannel: exponent B (in units of b) defined by A ~ Ae~ BW as 
a function of t'/t. Right pannels: Luttinger parameters K s 
(up) I<c (down) as function of ribbon width. 



reduced overlap diminishes the back-scattering terms and 
decreases the gap. Thus, larger values of t' render gap- 
less chiral-spin filtered states even in the presence of 
electron-electron interactions. In this regime, the sys- 
tem is described by models with no inter-edge Coulomb 
interactions [Tol 21 1. The data also show: K c < 1 and 
K s > 1 for all values of t' considered, indicating a strong 
tendency to form magnetically ordered structures at the 
edges. This conclusion is in agreement with ab-initio cal- 
culations that revealed strong magnetic instabilities (of 
ferro or anti- ferromagnetic nature) for ZGNR (23|. As 
Fig. (|3]) shows, a comparison with AGNR reveals gaps of 
an order of magnitude smaller in zigzag ribbons. Fur- 
thermore, the gap dependence A ~ 1/W (see Fig.(|4])) 
for AGNR changes to A ~ e~ w ; and reflects the dif- 
ferent nature of the respective edge states. For AGNR, 
the edge states are surface states (not localized along 
the edges) formed by a linear superposition of extended 
states near Dirac points [lij]. The I-SO interaction mildly 
localizes them along the edges with a localization length 
inversely proportional to the I-SO coupling [l^, 01 and 
makes them spin-filtered states. In contrast, the edge 
states in ZGNR are strongly localized in the absence of 
t' and the main effect of the I-SO interaction is to make 
them spin-filtered states with a suppressed unscreened 
Coulomb repulsion. 

In summary, we have obtained exact expressions for 
band-structure and wavefunctions of a ZGNR in its 
QSH phase, by solving a tight-binding Hamiltonian 
that includes intrinsic-spin orbit interactions. For ap- 
propriate ribbon widths, the I-SO interaction induces 
spin-filtered edge states with linear dispersion around 
k x = it. These states are composed by localized and 
damped-oscillatory wavefunctions and closely resemble 
those found in AGNR. Unscreened Coulomb interactions 



open a charge-gap in the spectrum while the spin-sector 
remains gaplcss as in AGNRs. However, numerical re- 
sults suggest an exponentially vanishing gap as a function 
of the ribbon's width W with gap amplitudes and decay- 
lengths decreasing with increasing values of the I-SO cou- 
pling t' . The exponential dependence of the gap is the 
first one reported for finite-width zigzag ribbons. Our es- 
timates give smaller values for gap amplitudes than those 
predicted in DFT calculations [l3( , while in general agree- 
ment with the trends observed in them. Furthermore, 
they emphasize the relevance of inter-edge unscreened 
Coulomb interactions in the QSH phase, as compared 
with previous results based on renormalization group ar- 
guments. The results also suggest an enhanced magnetic 
order along the edges reminiscent to the one found in 
AGNR for large values of the I-SO coupling [ljj. These 
results have immediate experimental consequences and 
could be tested in current setups 0]. 
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